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^ ! Abstract 

In the framework of the classical Glauber approach the exact analytical expression 
for the variance of the number of participants (wounded nucleons) for given centrality 
■ AA interactions is presented. It's shown, that in the case of nucleus- nucleus collisions 
rj^ ■ along with the optical approximation term the additional "contact" term appears. 

The numerical calculations for PbPb collisions at SPS energies show that at interme- 
diate values of the impact parameter the "optical" and "contact" terms contributions 
to the variance of the participants number are of the same order and their sum is in 
a good agreement with the results of independent MC simulations of this process. 
The correlation between the numbers of participants in colliding nuclei is taken into 
account. In particular it's demonstrated that in PbPb collisions at SPS energies 
the variance of the total number of participants approximately three times exceeds 
the Poisson one in the impact parameter region 10-12 Fm. The fluctuations of the 
O ' number of collisions are also discussed. 

a: 

jg 1 : 1 Variance of participants number in one nucleus 

At first we consider the variance V[iV,^(&)] of the number of participants (wounded nu- 
cleons) in one of the colliding nuclei N A (b) at a fixed value of the impact parameter b. 
In the framework of pure classical, probabilistic approach to nucleus-nucleus collisions, 
formulated in [I], we find for the mean value and for the variance of N£(b): 

(N£(b)) = AP(b) , (1) 

V[N£{b)] = AP(b)Q(b) + A(A - lW 12 \b) - Q 2 (b)} , (2) 
where P(b) = 1 - Q(b). For Q(b) and Q {12) (b) we have: 

Q(b) = J d ai T A { ai )[l - a^)] 3 , 
Q {12) (b) = J da 1 da 2 T A (a 1 )T A (a 2 )[l -ffi(ai) -0i(a 2 ) + a^ 12 \a u a 2 )] B , (3) 

where 

ci(ai,2) = / db 1 T B (b 1 )a(ai j2 - h + b) , 

<7< 12 >(ai, a 2 ) = j d6 1 T B (6 1 )(7(a 1 -b x + b)a(a 2 -h + b) . (4) 

Here T A and Tg are the profile functions of the colliding nuclei A and B ; the a (a) is 
the probability of inelastic interaction of two nucleons at the impact parameter value a 
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(/ a(a)da = <J™ N = a) and all integrations imply the integration over two-dimensional 
vectors in the impact parameter plane. 

The formula ([T]) and the first term in ([2]) correspond to the naive picture (so-called 
optical approximation) implying that in the case of AS-collision at the impact parameter 
b one can use the binomial distribution (|40p for N^(b) with some averaged probability 
P{b) of inelastic interaction of a nucleon of the nucleus A with nucleons of the nucleus B. 
At that the P(b) is considered the same for all nucleons of the nucleus A. 

The whole expression (J2J) is the result of more accurate calculation (see appendix lAl . 
when one uses probabilistic considerations taking into account the impact parameter plane 
positions of nucleons in the nuclei A and B and averaging then over these positions: 

V[N£(b)} = (N£(bf) - (N*(b)) 2 , (5) 

where 

. B A 

(X) = {{X) B ) A = XWTBib^dhWTA^dai . (6) 

k=l i=l 

Here X means average of some variate X at fixed positions of all nucleons in A and B; 
( ) A and ( ) B mean averaging over positions of these nucleons. 
In the limit <C Ra, Rb formulae (BJ reduce to 

0-1(01,2) ~ crT B (a lj2 + b) , 

<j( l2 \ai, a 2 ) ~ I(a 2 — cti) ■ T B (ai + b) , with 1(a) = j ' ds a(s) a(s + a) . (7) 

Note that in this limit the Q(b) and hence the first term of ([2]) and ([I]) depend only on the 
integral inelastic iViV cross-section <j = cr^^, but the Q^^(b) entering the second term of 
(J2J) depends also on the shape of the function a(b) through the integral 1(a) (see equation 

©). 

Note also that using of the approximation a (b) = a 5(b) for iViV interaction gives the 
same result (as taking of the limit <C Ra, Rb) only for naive part of the answer, which 
is expressed through Q(b). If one will try to use the approximation a(b) = a 8(b) for to 
calculate Q^ l2 \b), then one gets 1(a) = a 2 8(a) and = a 2 8(a 2 — ax)-Tg (ai + b) , which 
leads to infinite Q^ l2 \b) at B > 2. Meanwhile, for any correct approximation of cr(b), 
when <j(b) < 1 in correspondence with its probabilistic interpretation, we find definite 
finite answer for Q^ 12 \b). 

In our numerical calculations we'll use for NN interaction the "black disk" approxi- 
mation: 

<r(&) = 0(r*-|&|) , (8) 

or Gauss approximation: 

a(b) = exp(-& 2 /ri) . (9) 

In both cases a = nr%. For the nuclear profile functions Ta and Tb we'll use the standard 
Woods-Saxon approximation. 

We would like to emphasize that nontrivial second term in (j2J) arises only in the case 
of nucleus-nucleus collisions. For A = 1 or B = 1 it's equal to zero. At A = 1 due to 
explicit factor A — 1 in (j5J) and at B = 1 due to fact that in this case Q^ 12 \b) = Q 2 (b). 
This corresponds to the well known fact that for nucleus- nucleus collisions the Glauber 
approach doesn't reduce to the so-called optical approximation even in the limit 
Ra,Rb (see, for example, [2]). 
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PbPb, o=Jo(b)d 2 b=31mb, Ab=0 
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Figure 1: The variance of the number of wounded nucleons in one nucleus for PbPb colli- 
sions at SPS energies (cr = 0-^=31 mb, r N =l Fm) as a function of the impact parameter 
b. The points • and ■ - results of numerical calculations by formulae ([2]), ([3]) and ([7]) using 
respectively the "black disk" ([8]) and Gaussian (Q approximations for NN interaction; O 
and □ - results of independent MC calculations using for iViV interaction "black disk" ([8]) 
or Gaussian ([9]) approximation; * - "optical" approximation (the first term in formula 
([2])); + - the Poisson variance: V[N^{b)} = (N^(b)). The curves are shown to guide eyes. 



This additional term, which arises in ([2]) in the case of nucleus-nucleus collisions, 
depends, as we have mentioned, not only on integral value of inelastic NN cross-section 
a™ N = cr = / a (a) da, but also on the shape of the function a{b), i.e. on the details of iViV 
interaction at distances, which are much smaller than the typical nuclear distances. In 
quantum Glauber approach this corresponds to the fact that in the case of AA collisions, 
in contrast with pA collisions, the loop diagrams of the type shown in FigfS] appear. The 
typical momentum corresponding this loop integration is much larger than the typical 
nuclear momenta (which corresponds to smaller than typical nuclear distances) and one 
encounters the "contact" terms problem (see, for example, [21 El H]). The second term in 
formula (J2J) is the manifestation of this problem at the classical level. 
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Figure 2: An example of the loop diagram in quantum Glauber approach to AB-collisions. 
1 and 2 - nucleons of the nucleus A; 1' and 2' - nucleons of the nucleus B (see [21 El S] 
for details). 
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PbPb, o=Jo(b)d 2 b=31mb, Ab=0 
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Figure 3: The same as in Fig{IJ but for the mean number of wounded nucleons in one 
nucleus, calculated by formulae p]), d2J) and (j4j); * - "optical" approximation, calculated 
using formulae ([I]), ([3]) and (J7J). 

The numerical evaluation of the contribution of the additional - "contact" term in (T5]) 
are presented in FigfUfor PbPb collisions at SPS energies (er = (7^=31 mb, rjv=lFm). 
For the control we also carried out independent Monte-Carlo calculations of the mean val- 
ues and the variances involved presenting the results on the same figures. All calculations 
were done at fixed values of the impact parameter b (Ab = 0). 

In Fig{T]we see that the calculated "contact" term in ([2]) is essential and gives approx- 
imately the same contribution to the N^(b) variance for PbPb collisions at intermediate 
values of b, as the first "optical" term in (T5]). We see also that the results of indepen- 
dent MC calculations of the N^(b) variance are in a good agreement with the results of 
analytical calculations by formula ([2]), but only with taking into account its second term. 

Note that due to this "contact" term the N^{b) variance is larger than the Poisson 
one for peripheral PbPb collisions (at b > 7Fm). The week dependence of the results on 
details of NN interaction at nucleon distances is also seen. The results lay systematically 
slightly higher in the case of using the "black disk" ([8]) approximation for cr(b), than in 
the case of using the Gaussian (Q approximation with the same value of a . 

For the mean value {N^(b)), in contrast to the variance of N^(b), the exact answer 
coincides with the "optical" approximation result (see formula (JTJ and appendix and 
depends only on a. In Figj3]we see that MC calculations also confirm this result. 

2 Variance of the total number of participants 

Now we pass to the calculation of the variance of the total number of participants 
VWw(fy + ^(i)] at a fixed value of the impact parameter b. Clear, that we simply 
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PbPb, a=lo(b)d 2 b=31mb, Ab=0 
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Figure 4: The correlator between the numbers of wounded nucleons in colliding nuclei, 
calculated by formulae (fT2l) - (fT4l) and by independent MC simulations. The notations are 
the same as in Fig{TJ 



(10) 



have for the mean value 

{Ni(b) + N£(b)) = (N£(b)) + (NS(b)) 
and by (jHJ) for the variance 

+ ^(6)] = V[J^(ft)] + m« (&)] + 2{rt)^(6)) - (i^(6))(i^(6))} . (11) 
In naive approach ("optical approximation") there are no correlations: 

(Ndm5(b))-(N£{b))(N2{b)) = a. 

More accurate calculations (see appendix [B]) , based on formulae (jSJ) and (jB|), lead to 

(N£(b)NS(b)) ~ (N*(b))(N*(b)) = AB[QM(b) - Q(b)Q(b)} , (12) 

where 

Q{b) = ( da x T A { ai )[l - a x {ai)\ B , Q(b) = f ^^(61) [1 - ? 1 (& 1 )] A , 



QW(b) = J daATMaOTM&Otl - viiai)]*-^ ~ Mh)]^ ~ <?{ai ~h + &)] • (13) 
Here 

o-i(ai) = J db 1 T B (b 1 )a(a 1 - b\ + b) « aT B {ai + b) , 

?i(6i) = J da l T A {a 1 )a{a l -h + b)n aT A (b - b x ) (14) 
Note that Q(b) and o"i(ai) are the same as in formulae ([3]), (jlj) and (JTj). 
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PbPb, o=Jo(b)d 2 b=31mb, Ab=0 
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Figure 5: The same as in FigJH but for the variance of the total number of wounded 
nucleons N w (b) = N^(b) + N^{b) in colliding nuclei; the variance V[N w (b)} is calculated 
by formulae (T5]), ([3]), ([7]) and (|TTl) - (fl4T) . with taking into account the contribution of 



the correlator (N^(b)N^(b)) - (N£(b)){N%(b)) (see Fig0|; + - the Poisson variance: 
V[N w {b)] = (N w {b)}. 



The results of numerical calculations of the correlator ( fl2l) by formulae ( 1T31) and ( fl4|) 
together with the results of independent Monte-Carlo calculations for PbPb collisions at 
SPS energies are presented in FigJH 

Comparing Figj4] with Fig{l] we see that the contribution of this correlator to the 
variance of the total number of participants at intermediate values of b is about half of 
the variance for one nucleus V[iV^(&)] and is about the contribution of first "optical" term 
in (j2J). (Note that the relative contribution of the correlator to V[N^{b) + N^(b)] even 
greater at large values of b, b > 10.) The results are again in a good agreement with the 
results of MC calculations. 

In FigJH]we present the final results for the variance of the total number of participants 
in PbPb collisions at SPS energies, taking into account the contribution of this correlator. 
We see in particular that now the calculated variance of the total number of participants 
is approximately three times larger than the Poisson one in the impact 
parameter region 10-12 Fm. 

3 Discussion and conclusions 

In the framework of the classical Glauber approach the exact analytical expression for 
the variance of the number of participants (wounded nucleons) in AA collisions at a 
fixed value of the impact parameter is presented. It's shown, that along with the optical 
approximation contribution (which depends only on the total NN cross-section) in the 
case of nucleus-nucleus collisions there is the additional "contact" term contribution, 
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Figure 6: The same as in FigJSl but for the normalized variance V[N w (b)}/ (N w (b)) of the 
total number of wounded nucleons in colliding nuclei, N w (b) = N£(b) + N^(b). 



depending on the integral of the overlap of two inelastic NN cross-sections in the impact 
parameter plane. 

In the classical Glauber approach under consideration this "contact" contribution 
arises at correct taking into account the interactions between two pairs of nucleons in 
colliding nuclei (a pair in one nucleus with a pair in another). It's found, that the in- 
teractions of higher order, than between two pairs of nucleons, don't contribute to the 
variance. At that the expression for the mean number of participants proved to be exact 
already in optical approximation, based on taking into account only the averaged interac- 
tion between a nucleon in one nucleus with a nucleon in another. The same is also valid 
for the mean value and the variance of the number of NN-collisions in AA-interactions 
(see appendix |C]) . 

These results are obtained in the framework of the pure classical (probabilistic) Glauber 
approach [TJ. However it's possible to suppose, that in the quantum case the one loop 
expression for the variance and the "tree" expression for the mean number of participants 
and NN-collisions will be exact. 

Using the formulae obtained, the numerical calculations of the variance of the par- 
ticipants number in PbPb collisions at SPS energies are done. It's demonstrated that 
at intermediate values of the impact parameter the "optical" and "contact" term contri- 
butions are of the same order and their sum is in a good agreement with the results of 
independent MC simulations of this process. 

When calculating the variance of the total (in both nuclei) number of participants 
the correlation between the numbers of participants in the colliding nuclei is taking into 
account. The exact analytical expression for the correlator at a fixed value of the impact 
parameter is obtained. The results of numerical calculations of the correlator for PbPb 
collisions at SPS energies show that at intermediate and large values of the impact param- 
eter its contribution to the variance of the total number of participants is about half of 
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the variance in one nucleus, again in good agreement with independent MC simulations. 

In particular as a result it's found that the calculated variance of the total number 
of participants in PbPb collisions at SPS energies approximately three times larger than 
the Poisson one in the impact parameter region 10-12 Fm. (See FigsISl and [TUl for the 
normalized variance of the number of wounded nucleons and NN-collisions.) 

Note that the good agreement of the analytical and MC calculations ensures the relia- 
bility both of them and enables to use the developed MC algorithm in future experimental 
setup motivated calculations. 

The author thanks M.A. Braun and G.A. Feofilov for useful discussions. The work 
was supported by the grant RNP 2.2.2.2.1547 of Education Ministry of Russia and by the 
RFFI grant 06-02-16115a. 



Appendixes 

A Calculation of the participants variance for one 
nucleus 

The geometry of A5-collision is shown in FigJTl All a and b are the two-dimensional 
vectors in the impact parameter plane. In the framework of the classical (probabilistic) 
approach [1] the dimensionless a(b) is the interaction probability of two nucleons at the 
impact parameter value b: 

r a(b)db = a% l N = a (15) 



(all integrations imply the integration over two-dimensional vectors in the impact pa- 
rameter plane). This probabilistic interpretation means that a(b) < 1, so even in the 
limit rjy <C Ra,Rb we can't use approximation: a(b) = a5(b). The examples of valid 
approximations follows: the "black disk" approximation: 

a{b) = 6{r N - |6|) , a = irr 2 N , (16) 

the "grey disk" approximation: 

a{b) = j6{r N - \b\) , 7 < 1 , a = ^r 2 N , (17) 

Gauss approximation: 

_ b 2 

a(b) = Ce , C < 1 (!) , a = Cnr 2 N . (18) 

Ta and Tb are the profile functions of the colliding nuclei A and B. We'll imply the 
factorization takes place: 

A 

T A (ai, a A ) = II T A( a i) > ( 19 ) 

8=1 

which is sound approximation for heavy nuclei. Let us introduce also some shorthand 
notations: 

da^TA^da^ Jda^jTA^da^l. (20) 
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Figure 7: AB-collision. 



. B A , B A 

(X) = ((X) B ) A = XY[db k Y[d ai = / Xnr B (6 fe )d6 fc n r A(ai)dai (21) 

fc=i i=i ^ fc=i i=i 

= (X 2 ) - (X) 2 (22) 

Here X means average of some variate X at fixed positions of all nucleons in A and B; 
( ) A and ( ) B mean averaging over positions of these nucleons. 

We introduce now the set of variates X\, ...,Xa (each can take on a value equal only 
to or 1) by the following way: 

Xj = 1 if j-th nucleon of the nucleus A interacts with some nucleons of the nucleus B 

Xj = if j-th nucleon of the nucleus A doesn't interact with any nucleons of the nucleus B 

Then the number of participants (wounded nucleons) in the nucleus A in the given event 
can be simply expressed through these variates: 

tt = E^. (23) 

3=1 

So we have for the mean value: 



A A 

(« = EW = BW ( 24 ) 

i=i i=i 

and for the variance of N^{b): 

V[N£(b)] = (Ntibf) - (N*(b)) 2 , (N^bf) = (QT X,f) . (25) 

3=1 

Let us start our calculations from (^^ . Clear that for given configuration {a{\ and {bk}- 
X~ = 0- P{Xj = 0) + 1 • P(Xj = l)= Pj = l- Qj , (26) 
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where 

P(X j = 0) = q j =Y[(l-<r jk ) , (27) 
fe=i 

B 

P(X i = l)=p i = l-g J = l-n(l-^) , (28) 

fc=i 

ajk = cr(aj -bk + b) . (29) 

Here P{Xj = 0(1) ) is the probability that the variate Xj will be equal to or 1 corre- 
spondingly. We have to keep in mind that pj and qj are the functions of dj, bi,...,b B and 
b: 

Qj = <lj( a ji { & fc}, 6) , = Pj(aj, {h}, b) . (30) 

Straighforward calculations give: 

AAA A 

(N£{b)) = J2i(X]) B ) A = Y,«Pj)b)a = E((! - <&>b>a = A - E«<?;>b>a 
i=i i=i i=i j"=i 

B 

77 _ , TT / 77 _ , TT / 77. t . aJ 



(<?j)b = (LK 1 - a jfe))B = / II - = II / ^toC 1 _ = IK 1 ~ / ^fc^i 
fc=i fe=i fc=i fc=i 

= (1 - a,)* , a s = jdha n = J db^b^ -b 1 + b) = afa) (31) 



A 

((<1i)b)a = (^-^) B )a = /n^( 1 -°-i) fl = Jda j (l-a j ) B = j d aj T A { aj ){l- 
= Q{b) , = y rfaiT 4 (ai)(l-cr 1 ) B , ° x = \ dbiT B (b 1 )a(a 1 -b 1 +b) = ^(ai) 

A A 

(J\tf(6)> = A - Y,((h)b)a = A - E W = ^ - = A(l - Q(6)) = AP(fc) , (32) 

i=i i=i 

which coincides with formula ([T]) of the text. We see that the result for the mean value 
of the number of participants is the same as in an optical approximation. 
Let us now calculate the variance of N£(b). We start from ([25]): 

(N*(bf) = (ctxrf) = < £ x n x n ) = {: + e<4> ■ ( 33 ) 

i=i iu'2=i j'i^j'a=l i=i 

So we have to calculate the following two sums: 

E<*i> = E((^?)b)a = E<PQb>a (34) 

and 

A A A 

E (X h Xj 2 ) — ((Xj 1 X j2 ) B ) A = {{ x h ' Xj 2 ) B ) A . (35) 

A 

Note that the last expression can't be reduced to 7^ E ((^j'i)b)a ' ((-^72)5) a (0 

.71 7^2=1 

Just in this point the optical approximation breaks. 
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Clear that 

Xf = 2 • P{X j = 0) + l 2 • P(X j = 1) = X] = Pj = 1 - q 3 . 
Notations are the same as in (1271) - (|32j) . And for the first sum (1341) we find: 

A A A 

E<*i> = E«*S%>a = J2(C^)b)a = (N£(b)) = AP(b) = A(l - Q{b)) . (36) 
j=i j=i j=i 

For the second sum (I3^|) the straighforward calculations give: 



x h x h = x h ■ x h = VhVh = 0- - - Qk) = 1 - 9ii - fca + ^1*2 > 

E ( X ii X i 2 ) = E (( X ~ fci ~ ^2 + Vh^BlA = 
ii^i2=i ii^i2=i 

(yl A \ A 

E<foi>B>A+£«<&>B>A + E 
Jl=l J2=l / jl^j2 = l 

= A(A-i)-(A-i)A(g(6) + g(6))+ E ({QjiQfa) b) a ~ 

31^32=1 

= A(A-l)[l-2Q(b) + Q^\b)] , 

where we introduce 

^ ' ii^j2=i 

Let us now calculate Q^ 12 ^(6): 



(ihQj^B = (IK 1 - a hhi) IK 1 ~ °i 3 *a))B = / II - OjifcX 1 - a i2fc) = 
fci=i fc 2 =i fc=i 

y d&i(l - cr^i - a j2 i + cr^ia^i)^ = (l - J db\o hl - J db x o h x + J dbxa hl a 32 i 



= (1 _ aji - aj2 + a ^) B , 
where o 3x and a 32 are given by (I3TI) and 

a (jih) = J db\<jj x i<jj 2 \ = J dbiT B {bi)a{aj 1 — &i + b)a{a 32 — b\ + b)= a^ l2 \aj x ,aj 2 ) (37) 
So for Q( 12 )(6) we find 

Q (12) (^ I r^r TT E «^ 2 >b>a = jr^s £ ((i - ^ - ^ + ^ li2) ) B >A 

^rry . e / n ^(i - - + 



a(a-i) 



Ji^i2=l »=1 
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1 A r * 

= A(A-l) ^ / da n da ^ 1 ~ a n ~ a h + a° lj2) ) B 

= J da 1 da 2 T A (a 1 )T A (a 2 )(l - <Tl -a 2 + a (12) ) B , (38) 

where 

cr (12) = / dbi(Jn(J2i = J dbiT B (bi)a(ai - bi + b)a(a 2 - h + b) = a ( - 12 \ai,a 2 ) . 
Substituting ([3"3l-(l3"8l into ([25]) we find for the variance of N£(b): 

A A 

V{N*(b)} = (N^bf) - <iV^)> 2 = £ (%-X*> + EW - <^( fe )> 2 

j"i^j2=i i=i 

= A(A - 1)[1 - 2Q(6) + Q (12) (^)] + A(l - Q(6)) - [A(l - Q{b))} 2 

= AQ(b) - A 2 Q 2 (b) + A(A - \)Q^ 2 \b) = AQ(b)[l - Q(b)] + A(A - l)[Q {12) (b) - Q 2 (b)] , 

(39) 

which coincides with the formula (j2j) of the text. 

The naive approach (optical approximation) implies the binomial distribution for N A : 

P opt {N*) = Cf P(bf™ Q(b) A - N ", P(b) = 1 - Q(b) (40) 
which immediately leads to 

(N*(b)) opt = AP{b) (41) 

and 

V[N^{b)U = AP{b)Q{b) = (Ni{b))Q{b) . (42) 

We see that this gives true answer only for the mean value (N^(b)). For the variance 
V[N^(b)} the results coincide only at A = 1 due to explicit factor A — 1 in (1351 or at 
B = 1 as in this case Q^ l2 '{b) = Q 2 {b) (i.e. for pA-collisions). For nucleus- nucleus 
collisions, when A > 2 and B > 2 the naive result (142]) for variance is not valid (see text 
for details). 

Note that for peripheral AA collisions (at large b), when P(b) becomes small (P(b) <Cl, 
Q(b) ~ 1), the naive distribution (HOI) and the variance (l4"2"j) reduce to the Poisson ones: 
VWMU = W(6)> (^e FigflJ. 

B Correlation between the numbers of wounded nu- 
cleons in colliding nuclei at fixed centrality 

The calculations are similar to ones in the appendix [A] (we use the same notations). 
Allong with the set of variates X±,...,X A we introduce in the symmetric way also the 
set of variates X\, ...,Xb (each can again take on a value equal only to or 1) by the 
following way: 

Xk = 1 if k-th nucleon of the nucleus B interacts with some nucleons of the nucleus A 
Xk = if k-th nucleon of the nucleus B doesn't interact with any nucleons of the nucleus A 
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Then similarly to (1231 the number of participants (wounded nucleons) in the nucleus B 
in the given event can be simply expressed through these variates: 

N B w {b) = jzX k . (43) 
k=i 

Then 

(N£(b)NM) = E EPW = J2J2((x~Z) B ) A m 

3=1 k=l j=l k=l 

and 

XjX k = 1 • 1 • P jk {l, 1) + 1 • • 0) + • 1 • P jk (0, 1) + • • P iJfe (0, 0) = P jk {l, 1) . (45) 

Here Pj k (0(l), 0(1)) is the probability that the variates Xj and X k will be equal to or 1 
correspondingly. We have 

Pjfe(l, 1) = CTjfe + (1 - cr jk )p jk p jk , (46) 

where cr^ is the probability of the interaction of the j-th nucleon of the nucleus A with 
the fc-th nucleon of the nucleus B: 

a jk = cr(aj - b k + b) , (47) 

and pjk is the probability of the interaction of the j-th nucleon of the nucleus A with 
at least one nucleon of the nucleus B except the k-th nucleon (correspondingly pj k is 
the probability of the interaction of the fc-th nucleon of the nucleus B with at least one 
nucleon of the nucleus A except the j-th nucleon): 

B A 

= i - n ( x ~ a jk') , pjk = i - n i 1 ~ a j'k) ■ ( 4§ ) 

k'=i(k'^k) j'=i(j'¥=j) 

Combining (Jli|) - (j4^]l and acting as in the appendix |X] we find the formulae (fT2]) - (fT4"|) 
of the text. 



C On fluctuations of the number of collisions 

In this appendix we discuss briefly the fluctuations of the number of NN-collisions in 
AA-interactions at fixed value of centrality in the framework of the approach under con- 
sideration. 

To calculate the number of collisions we define the set of A variates Y\_ , . . ., Ya (each 
can take on a value from to B) by the following way: 

Yj : = if j-th nucleon of the nucleus A doesn't interact with any nucleons of the nucleus B 
Yj = 1 if j-th nucleon of the nucleus A interacts with one nucleon of the nucleus B 
Yj = 2 if j-th nucleon of the nucleus A interacts with two nucleons of the nucleus B 
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Figure 8: The mean number of NN-collisions in PbPb interactions at SPS energies calcu- 
lated by formulae (1551) and (I56I) as a function of the impact parameter b; * - "optical" 
approximation, calculated using formulae (1571) and (|59|) . The notations are the same as 
in FigHJ 



Yj = B if j-ih nucleon of the nucleus A interacts with all nucleons of the nucleus B 

Then N c (b) (the number of NN-collisions in the given event with impact parameter b) can 
be expressed through these variates as follows: 

N c (b) = J2Y J (49) 

Clear that again (see appendix fA|) : 

P{Y = 0) = P(X j = 0) = gj = n (1 - oi-fc) (50) 

k=l 

To calculate P(ij = n) forn = 1, B let us introduce {&i, A;„} - the random sampling 
from the set {1, B} and {k n +i, ks} - the rest after sampling. Then 

P(Yj = n)= Yl Vjki-Vjkni 1 ~ - Oiks) (51) 

{fci,...,fe n } 

We can calculate the mean value of the number of collisions: 

A A 

(N c (b)) = £<Y}> = £«^)b)a (52) 
i=l i=l 

For the given configuration {aj} and {bk} we have: 

^=EnP(K i = n) (53) 

n=0 
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Figure 9: The variance of the number of NN-collisions in PbPb interactions at SPS energies 
as a function of the impact parameter b; * - "optical" approximation, calculated using 
formulae (1571) and (|60|) ; + - the Poisson variance: V[N c (b)} = (N c (b)). The notations 
are the same as in Fig{TJ 



and 



B B 

n=0 {ki,...,k„} 



n=0 



B 



(J; 



n=0 " fe=l {ki,...,k n } 
B 



= ^C; ff »(l- ffj f» = Ba,. 

n=0 

We use the same notations (see (IHTj) ) as in appendix [XI 

Oj = y db\Oj\ = J dbiT B {bi)a{aj — b\ + 6) ~ aT B (aj + 6) . 
Finally we find: 

/I yl A A A A 

j = l j = l j = l j = l j=l 8=1 



(54) 



(55) 



where 



= / daifTi = / cfaidfoifTn = / daid&iT^aijT^&iJa^ai — &i + b) (56) 
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Figure 10: The same as in FiglHl but for the normalized variance, V[N c (b)}/ (N c (b)) , of 
the number of NN-collisions. 



and at r N <C Ra, Rb 

X(b) ps cr J da 1 T A (a 1 )T B (a 1 + b) , (57) 

which coincides with the optical approximation result. 

Really, assuming the binomial distribution for N c (b) with the averaged probability 
x{b) of NN-interaction: 

P op t(N c ) = Cf B X (b) N ° [1 - X (b)] AB - Nc , (58) 

we have 

(N c (b)) opt = AB X (b) (59) 

and 

V[N c (b)] opt = AB X (b)[l - X (b)) = (N c (b))[l - x(b)) . (60) 

Note that for heavy nuclei x{b) is small even for central collisions {x(b) ~ r 2 N jR\ ^Cl), 
and the distribution f)58p and the variance (IBT)]) practically coincide with the Poisson ones: 
V[N c (b)} opt ~ (JV C (6)>. 

Comparing fl55|) and fl59l) we see that the optical approximation gives true answer for 
the mean value (N c (b)). Although we could not find the closed formula for the variance 
V[N c (b)], which would be similar to formulae ([2]), ffTTT) and f[T2]) for the variance of the 
number of participants V[N^(b)} and V[N^{b) + (b)], our calculations show that the 
expression for the variance of N c (b) will again differ from the optical approximation one 
(|60|) . The results of our separate Monte-Carlo simulations confirm these conclusions. 

In Figs JH|9l and [TU1 we present the results of our calculations of the mean value (N c (b)) 
and the variance ^[^(6)] in PbPb collisions at SPS energies at different values of the 
impact parameter b. 
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